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ABSTRACT 



Understanding neural network dysfunction in neurodegenerative disease is imperative to effectively develop 
network-modulating therapies. In Alzheimer's disease (AD), cognitive decline associates with deficits in 
resting-state functional connectivity of diffuse brain networks. The goal of the current study was to test whether 
specific cognitive impairments in AD spectrum correlate with reduced functional connectivity of distinct brain 
regions. We recorded resting-state functional connectivity of alpha-band activity in 27 patients with AD spec- 
trum — 22 patients with probable AD ( 5 logopenic variant primary progressive aphasia, 7 posterior cortical atro- 
phy, and 10 early-onset amnestic/dysexecutive AD) and 5 patients with mild cognitive impairment due to AD. 
We used magnetoencephalographic imaging (MEGI) to perform an unbiased search for regions where patterns 
of functional connectivity correlated with disease severity and cognitive performance. Functional connectivity 
measured the strength of coherence between a given region and the rest of the brain. Decreased neural connec- 
tivity of multiple brain regions including the right posterior perisylvian region and left middle frontal cortex cor- 
related with a higher degree of disease severity. Deficits in executive control and episodic memory correlated 
with reduced functional connectivity of the left frontal cortex, whereas visuospatial impairments correlated 
with reduced functional connectivity of the left inferior parietal cortex. Our findings indicate that reductions in 
region-specific alpha-band resting-state functional connectivity are strongly correlated with, and might contrib- 
ute to, specific cognitive deficits in AD spectrum. In the future, MEGI functional connectivity could be an impor- 
tant biomarker to map and follow defective networks in the early stages of AD. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creafivecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Synchronization of neuronal activity across different brain regions is 
a fundamental property of cortical networks (Singer, 1999). Spontane- 
ous correlations in activity occur between brain regions that are spatial- 
ly distributed yet functionally related (Fox and Raichle, 2007). Patterns 

Abbreviations: CDR-SOB. Clinical Dementia Rating Sum of Boxes: CVLT, California 
Verbal Learning Test: fMRl, functional magnetic resonance imaging: IvPPA. logopenic 
variant primary progressive aphasia: MCI, mild cognitive impairment: MEGI, 
magnetoencephalographic imaging: MMSE, Mini-Mental State Exam: PCA. posterior 
cortical atrophy: VOSP, Visual Object and Space Perception. 
•• Correspondence to: Professor in Residence, Department of Radiology and Biomedical 
Imaging, UCSF School of Medicine, University of California San Francisco, 513 Parnassus 
Avenue, S362, San Francisco, CA 94143. 

* Correspondence to: Assistant Professor, Gladstone Institute of Neurological Disease, 
UCSF Memory and Aging Center , 1 650 Ovirens St, San Francisco, CA 941 58. 

E-mail address: Srikantan.Nagarajan@ucsfedu (S.S. Nagarajan), 
KVossel@memoi'y.ucsfedu (KA. Vossel). 



of temporally coherent fluctuations within widely distributed cortical 
systems also correspond with behavior (Hampson et al., 2006; Seeley 
et al., 2007). Large-scale brain networks that are developmentally 
formed to regulate behavior are also selectively targeted by different 
neurodegenerative diseases (Seeley et al., 2009). 

Alzheimer's disease (AD), the most common age-related dementia, 
is associated with early disruption of functional neural networks. Struc- 
tural, molecular, and metabolic changes in the AD brain mirror, to a 
large extent, the distributed intrinsic network collectively called the 
default mode network — regions characterized by task-related deacti- 
vation on resting-state functional magnetic resonance imaging (fMRl) 
(Buckner et al., 2005; Seeley et al., 2009). Regions that show the stron- 
gest functional connectivity, such as the default mode network, are 
often the most vulnerable to neurodegenerative disease (Seeley 
et al., 2009; Zhou et al., 2012). Electroencephalographic (EEC) and 
magnetoencephalographic (MEG) recordings have demonstrated reduc- 
tions in global neural synchrony in AD and mild cognitive impairment 
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(MCI) (de Haan et al., 2008; Jelic et al., 1996; Stam et al., 2009). These 
changes in synchrony may indicate failing connectivity within large- 
scale functional networks. The relationship between region-specific net- 
work connectivity in AD and domain-specific cognitive dysfunction, how- 
ever, remains largely unknown. 

In this study, we tested the hypothesis that region-specific resting- 
state functional connectivity patterns are correlated with focal cognitive 
deficits in AD spectrum. We explored resting-state functional connec- 
tivity using MEG imaging (MEGI)-acquired alpha-band oscillations in 
patients meeting research criteria for probable AD or MCI due to AD. 
Our results demonstrate the differential involvement of distinct function- 
al networks and the heterogeneous nature of AD. Robust correlations to 
cognitive deficits indicate region-specific functional connectivity deficits 
as a potential tool to gauge specific network dysfunctions in AD. 

2. Methods 

2.1. Subjects 

Patients were recruited from research cohorts at the University of 
California San Francisco (UCSF) Memoiy and Aging Center, and consisted 
of 27 patients with AD spectrum, including five with logopenic variant 
primary progressive aphasia (IvPPA), seven with posterior cortical atro- 
phy (PCA), ten with amnestic/dysexecutive variant, and five with MCI 
(four amnestic MCI and one language-predominant MCI). All patients 
undei"went a complete clinical history and physical examination, a neuro- 
psychological assessment, and a structured caregiver interview. Diagnosis 
was made by consensus at a multidisciplinary meeting and was based on 
clinical evaluation, neuropsychological test performance, and biomarkers 
(Supplementary Table 1 ). Of the 23 patients with amyloid imaging and/or 
CSF analysis for beta-amyloidi^2 peptide and tau protein, 21 had findings 
supportive of AD. The two subjects that were amyloid imaging negative 
(one borderline) were both MCI. Each had hippocampal atrophy and def- 
icits in episodic memoiy, which are strong predictors of conversion to AD 
(Devanand et al., 2007). Forty-eight percent (11 of 23) of the AD- 
spectrum cohort with apolipoprotein E (apoE) genotyping were carriers 
of the apoE4 allele. Patients received the initial diagnosis of probable AD 
or MCI if they fulfilled the National Institute of Aging-Alzheimer's Associ- 
ation criteria for probable AD or MCI due to AD, respectively (Albert et al., 
2011; McKhann et al., 2011). Each probable AD presentation was then 
reviewed to assess whether patients fulfilled the specific diagnostic 
criteria for PCA or IvPPA (Gorno-Tempini et al., 2011; Mendez 
et al., 2002). Patients who met the criteria for probable AD but did 
not meet criteria for either PCA or IvPPA were categorized into 
amnestic/dysexecutive variant. We also recruited 15 healthy age-, 
education-, and sex-matched control subjects (mean age 64.3 ± 
5.0, p = 0.59, mean education 16.8 ± 2.2, p = 0.18, and mean 
Mini-Mental State Exam score 29.5 ± 0.7, p < 0.0001 vs. patients, un- 
paired t tests; percent male 53%, p = 0.43, Chi-square test) for the com- 
parisons of MEG power spectrum analysis and gray matter volume 
assessment based on structural magnetic resonance imaging (MRI). In- 
formed consent was obtained from all participants or their assigned 
surrogate decision makers. The study was approved by the UCSF Com- 
mittee on Human Research. 

2.2. Neuropsychological assessment and statistical analysis 

Each patient completed the Mini-Mental State Exam (MMSE), Clini- 
cal Dementia Rating scale (CDR) and CDR Sum of Boxes (CDR-SOB) 
(Morris, 1993) (Table 1), and a batteiy of neuropsychological tests de- 
signed to assess executive function, spatial ability, memory, and lan- 
guage function (Table 2). Tasks used to assess the individual domains 
of cognitive function were detailed in a previous report (Kramer et al., 
2003) and included D words, animals, digit span backward, a modified 
Trail Making Test, Stroop color naming, Stroop interference, Boston 
Naming Test (BNT), California Verbal Learning Test-Short Form 



Table 1 

Patient demographics. 



Age (years) 


63.2 ± 7.1 


MMSE (/30) 


22.5 ± 4.2 


Sex (% male) 


41 


Handedness (% rigiit) 


92 


Education (years) 


15.8 ±2.2 


CDR 


0.84 ± 0.3" 


CDR-SOB 


5.1 ± 1.8* 


Age at onset (years) 


57.8 ± 7.1 


Disease duration 


5.4 ± 3.3 



Tile patient group includes a total of 27 patients with posterior cortical atrophy (n = 7), 
amnestic/dysexecutive (n = 10), logopenic variant primaiy progressive aphasia (n = 5), 
and mild cognitive impairment (n = 5). Numbers show means and standard 
deviations. 



CDR = Clinical Dementia Rating. CDR-SOB = Clinical Dementia Rating Sum of Boxes, 
MMSE = Mini-Mental State Exam. 
" n = 26. 



(CVLT) containing 9 items, CVLT 30-second free recall, CVLT 10-minute 
free recall, Benson copy, Benson free recall, and the number location 
test of the Visual Object and Space Perception (VOSP) batteiy. We catego- 
rized the fluency tasks under executive function although we acknowl- 
edge that they also rely on language (Possin et al., 2014). Given that 
fewer than 25 patients performed the modified Trail Making Test and 
Stroop test, we excluded these tasks from the subsequent correlation 
analyses (Cohen, 1992). We performed one-way-ANOVA with Tukey 
post hoc tests to compare the cognitive scores between subgroups of AD 
spectrum using SAS (SAS Institute, Caiy, NC). 

2.3. Magnetic resonance image acquisition 

Structural images were acquired on a 3 Tesia Siemens MRI scanner at 
the Neuroscience Imaging Center-UCSF for 22 of the 27 patients and for 
all 15 control subjects. The remaining 5 patients were studied using 
MRI scans obtained at an outside facility within 2 years of their MEG 
evaluation. 

2.4. Magnetoencephalographic image acquisition and processing 

MEGI uses MEG sensor data with millisecond precision and applies 
source reconstruction algorithms to overlay cortical oscillatory activity 
onto structural brain images. This process enables reconstruction of 
the oscillatoiy activity from MEG data in specific brain regions with 
high spatiotemporal precision. 

Each subject unden/vent MEGI on a 275-channel whole-head MEG 
system (MISL, Coquitlam, British Columbia, Canada) consisting of 275 
axial gradiometers (sampling rate = 600 Hz). Three fiducial coils in- 
cluding the nasion and the left/right pre-auricular points were placed 
to localize the position of the head relative to the sensor array. These 
points were later co-registered to a Tl -weighted MRI to generate a 
head shape. Data collection was optimized to minimize within-session 
head movements to not exceed 0.5 cm. Ten minutes of continuous re- 
cording was collected from each subject lying supine, awake, and with 
eyes closed. We selected a 60-second epoch of contiguous stationary 
segment of data for MEG source data analysis. The full recording of the 
data was not used to minimize the inclusion of motion artifact and 
other potential artifacts. Previous studies have determined that a 60- 
second window provides reliable, consistent power for reconstruction 
of brain activity from the resting-state MEG data (Guggisberg et al., 
2008; Hinkley et al., 2010; Hinkley et al., 2011). Artifact free segments 
were selected by using signal amplitudes less than a threshold of 10 
pT, and by visually inspecting the sensor data to identify segments with- 
out artifacts generated by eyeblinks, saccades, head movements, or 
muscle contractions. Both controls' and patients' data were within our 
specified limits of signal scatter, and the two groups had no difference 
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Table 2 

Neuropsychological performance scores. 



Ncuropsychologicsl test 






PCA 


Arnn/Dys 


IvPPA 


MCI 


P 


CDR-SOB 




26 


5.6 ± 1.6 


5.3 ± 1.9 


4.6 ± 1.1 


4.4 ± 2.4 


0.54 


Executive function 


Lexical fluency 


27 


12.1±3.2 


9.2 ± 2.3 


6.2±3.2 


8.8 ± 2.3 


<0.05 




Category fluency 


27 


10.0 ±5.1 


11.3 ±4.2 


6.8 ± 3.6 


11.2 ± 5.3 


0.40 




Modified Trails (lines/s) 


22 


0.06 ± 0.0 


0.2 ± 0.2 


0.15 ± 0.1 


0.3 ± 01 


0.21 




Stroop color naming 


21 


26±16 


57 ±25 


34 ±15 


64±14 


<0.05 




Stroop inhibition 


21 


10 ± 8.1 


23 ± 15 


10 ± 8 


32 ± 13 


< 0.05 




Digits span backward 


27 


2.6±1.3 


3.1 ± 0.7 


2.5±0.6 


4.4±1.1 


< 0.05' 


Language 


Boston Naming Test 


27 


13 ± 5.5 


12 ± 3.0 


10 ±3.1 


11 ± 1.7 


0.68 


Visuospatial 


Benson copy 


27 


3.8±2.1 


13.1 ±4.0 


13.7±1.3 


15.2±0.8 


< 0.0001 " 




VOSP number location 


25 


4.1 ±1.0 


7.4±2.6 


7.2 ± 3.2 


9.4±0.5 


< 0.01" 


Learning & memory 


CVLT total score 


25 


17.2 ± 5.5 


17.4 ± 5.8 


17 ±5.3 


20 ± 5.3 


0.909 




CVLT 30 s recall 


25 


4.0 ± 2.7 


3.4 ± 1.1 


5 ±3.1 


4± 1.0 


0.569 




CVLT 10 m recall 


25 


2.8 ± 2.7 


1.2 ± 1.4 


4.2 ± 3.4 


2.2 ± 1.7 


0.158 




Benson 10m recall 


27 


1.8 ± 1.8 


1.8 ± 2.9 


5.8 ±3.2 


5.6 ± 4.6 


<ao5 



p values indicate the significance of one-way-ANOVA. 

Bold text indicates the groups that were statistically different after Tukey post hoc comparison between the patient groups of posterior cortical atrophy (PCA). logopenic variant primary 
progressive aphasia (IvPPA). amnestic/dysexecutive (Amn/Dys). and MCI. with the significance threshold set to 0.05. 

CDR-SOB = Clinical Dementia Rating Sum of Boxes, CVLT= California Verbal Learning Test, VOSP number location = number location task of Visual Object and Space Percepfion battery, 
30 s = 30 seconds, 10 m = 10 minutes. Shown are means and standard deviations. 
" MCI vs. IvPPA and PCA 

PCA vs. MCI, IvPPA, and Amn/Dys groups. 

PCA vs. MCI and Amn/Dys groups. 



in head movement (unpaired t test comparing tlie degree of tiead 
movement before and after the resting data segment. Mean ± SD 
for patients = 0,34 ± 0.28 cm, controls = 0.27 ± 0.18 cm; t = 0.85, 
p = 0.4). 

In the current study we used a 20-mm isotropic grid in the sensor 
space, and reconstracted into source space using a minimum-variance 
adaptive spatial filtering technique (Dalai et al., 2008). Tomographic re- 
constructions of the MEG data were generated using a head model 
based on each subject's structural MRI. Each subject's Tl -weighted 
MRI was spatially normalized to the standard Montreal Neurological 
Institute template brain using standardized procedures (5 mm voxels; 
SPM2 http://www.fil.ion.ucl.ac.uk/spm/software/spm2) (Ashburner 
and Friston, 1999). The results of normalization were manually verified 
in all subjects. A whole brain volume of interest for lead field computa- 
tion was generated by transformation of all the points within a spatially 
normalized MRI that corresponded to locations within the brain and ex- 
cluded non-cerebral points. Source-space MEGl analysis and functional 
connectivity were generated using the NUTMEG software suite ( http:// 
nutmeg.berkeley.edu) (Dalai et al., 2011). For source-space recon- 
structions, the MEG sensor data were filtered with a phase-preserving 
bandpass filter (fourth-order Butterworth; 1-20 Hz). The source-space 
reconstmction approach allowed amplitude estimates at each voxel de- 
rived through a linear combination of a spatial weighting matrix with 
the sensor data matrix (Hinkley et al., 201 1 ). The alpha frequency peak, 
identified as the greatest power density within the 8-12 Hz range for 
each subject was chosen to select a subject-specific alpha band. Mean 
alpha power of a 4-Hz window around this peak was used for the subse- 
quent functional connectivity analysis. Previous studies have shown 
that, out of a range of frequency bands, the alpha-band connectivity 
maps provide the greatest within-session and cross-session reliability 
measurements and highest signal-to-noise ratio (Guggisberg et al., 
2008; Hanslmayr et al., 2011; Hindriks et al., 2011; Hinkley et al., 
2012; Hinkley et al., 2011; Klimesch, 1999; Palva et al., 2010; 
Westlake et al., 2012). 

We computed imaginary coherence, which is a reliable metric for 
functional connectivity with MEG reconstruction (Engel et al., 2013; 
Guggisberg et al., 2008; Martino et al., 2011; Nolte et al., 2004). 
Coherence-measures quantify the strength of neural oscillations. An im- 
portant limitation of coherence is the introduction of signal mixing arti- 
facts, where a single active source contributes to the other sensors 
through volume spread. We chose imaginary coherence because it 
captures only the coherence that cannot be explained by volume 



spread, and thus gives more accurate results. Previous studies have 
demonstrated that imaginary coherence reduces overestimation biases 
of MEG data and is able to sample interactions between source time- 
series, independent of the class of spatial filter used (Guggisberg et al., 
2008; Martino et al., 201 1 ; Nolte et al., 2004). In addition, more recent 
studies utilizing MEG in preoperative assessments have demonstrated 
the biological significance of imaginary coherence in evaluating func- 
tional connectivity (Martino et al., 201 1 ; Tarapore et al., 2012). A single 
voxel's Fisher's Z-transformed imaginary coherence values between 
that voxel and all other voxels in the grid were averaged to compute 
the global connectivity at each voxel. Global connectivity, repre- 
sented by the mean imaginary coherence measure in our study, 
would be akin to mean path-length of the maps of graph theoretic 
measures. 

Neuropsychological measures from each patient, including the CDR- 
SOB and cognitive test scores, were correlated with global connectivity 
values at each voxel using the Pearson's correlation. For all voxel-wise 
functional connectivity correlations with neuropsychological scores 
that survived an uncorrected p < 0.05, we performed a multiple compar- 
isons correction using a False Discovery Rate (FDR) set at 10% 
(Benjamini and Hochberg, 1995). When effects were significant, we ex- 
plored how robust these effects were by applying a more stringent 
threshold, with an FDR cutoff value of 5%. We report the most conserva- 
tive multiple comparisons correction producing significant effects. Con- 
trol subjects were excluded from this analysis because their high scores 
and low variability on the neuropsychological measures created ceiling 
effects that precluded such correlations. Group comparison of spectral 
power between the patient group and the control group was done by 
using an unpaired t-test. 

2.5. Voxel based morphometry 

The voxel based morphometry (VBM) analysis used subjects with 3 
Tesla MRI scans (n = 22 patients and 15 controls) to maintain consis- 
tency across subjects and to generate quality data for atrophy pattern 
assessments. We assessed the gray matter volume variations of patients 
and controls using standard methods of optimized VBM procedures 
within the VBM8 toolbox of statistical parametric mapping version 8 
(SPM8) ( http://www.fil.ion.ucl.ac.uk/spm/software/spm8/ ). SPM8 is 
equipped with an Expectation Maximization Segmentation tool for 
tissue segmentation to obtain probabilistic maps of gray matter and 
white matter. After tissue segmentation, gray matter probabilistic maps 
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Fig. 1. Averaged power spectral density estimates for different frequency bands. AD-spectrum 
patients showed decreased power of higher frequency activity over the alpha (p < 0.01 ), 
beta (p < 0.01 ) and gamma range (p = 0.08), and an apparent increase of power in 
lower frequency activity over the delta and theta range (nonsignificant), compared to 
age-matched healthy controls. Alpha range is highlighted in yellow. Shaded zone around 
each line depicts standard error. Spectral data were derived from MEG sensors, n = 27 
patients, n = 15 controls. 



of all subjects were initially transformed into the standard Montreal 
Neurological Institute space. The segmented gray matter probabilis- 
tic maps in their respective native space were then spatially normal- 
ized again to the population-based templates using a nonlinear 
transformation. Jacobian modulation was performed to obtain the 
volume differences from the spatially normalized gray matter 
image. The modulated images were smoothed with a 12-mm full- 
width-at-half-maximum isotropic Gaussian kernel to be used for 
group analyses. 



Group comparisons between patient subgroups and the control 
group were done using unpaired t-tests. Significant group differences 
were identified at family- wise error p < 0.05 unless otherwise indicated. 
Age was included as a nuisance covariate. We also performed an analy- 
sis in which we correlated each of the voxel-wise volume estimates 
with the neuropsychological performance scores, while keeping the 
total intracranial volume and age as covariates. Significant clusters of 
correlation were determined using height and extent thresholds of 
p < 0.001 uncorrected, at cluster level. 

2.6. Bayesian validation analysis 

A Bayesian hierarchical model was used to further validate the re- 
sults from the correlation analysis. This analysis considered the signifi- 
cantly correlated voxels within an immediate cluster around the peak 
voxel (i.e., voxel that shows the strongest correlation). To maintain 
consistency of voxel-cluster selection across different cognitive tasks, 
we kept it limited to the immediate cluster of 27 connected voxels in 
three-dimensional space. For comparability across tests, we standard- 
ized the scores from each test (by subtracting off the mean score and di- 
viding by the standard deviation). Let yivt denote the standardized score 
of subject i corresponding to voxel v on test t. The likelihood of our 
model was given by: 

Vive = Pvt^iv + + ^ivt > 

where (ivt is the effect on the score of test t of the connectivity (imag- 
inary coherence) (Xi„) of voxel v. The random effect b, captures the cor- 
relation of repeated measurements on subject i. The error term 8ivt 
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(max) 



D r = -0.59, p < 0.01 

(max) 




10 
8 
6 
4 
2 



-If. 



• MCI 

■ AD-Amn/Dys 

^ AD-Language 

A AD-PCA 



0.04 0.06 0.08 

Imaginary Coherence 



0.04 0.06 0.08 

Imaginary Coherence 



Fig. 2. Reduced resting-state functional connectivity is correlated with disease severity in multiple brain regions in AD spectrum. Imaginary coherence, reflecting the resting-state alpha- 
band global functional connectivity, predicts the degree of disease severity, as measured by the clinical dementia rating sum of boxes (CDR-SOB). The voxel containing the peak correlation 
between CDR-SOB and imaginary coherence identified over the right posterior perisylvian region (indicated by the "X" in panel (A) and plotted in panel (C)),and over the left frontal cor- 
tex (indicated by the "X" in panel (B) and plotted in panel (D)) are shown. Voxel-wise multiple comparisons are thresholded with 5% FDR correction. Statistical maps are superimposed on 
a rendering of the Montreal Neurological Institute template brain. The color scheme is normalized to the peak voxel correlation, p values on scatter plots indicate corrected p value for the 
peak voxel coixelation. n = 27 patients. MCI = mild cognitive impairment, AD = Alzheimer's disease, PCA = posterior cortical ati ophy, Amn/Dys = amnestic/dysexecutive, AD-Language = 



logopenic variant primary progressive aphasia, CDR-SOB = 
the peak voxel. 



Clinical Dementia Rating Sum of Boxes, i^oae 



-- correlation coefficient at the 5% FDR threshold, r(„ 



: correlation coefficient of 
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Fig. 3. Resting-state functional connectivity deficits in the left dorsolateral prefrontal cortex correlate with impairments in cognitive performance in AD spectrum. Resting-state functional 
connectivity as measured by imaginaiy coherence of the left dorsolateral prefrontal cortex correlated with performance of (A, F) lexical fluency (D words), (B, G) category fluency (ani- 
mals), (C, H) digit span bacl<ward. (D, 1) CVLT 30-second recall, and (E, J) CVLT total score. Performance on CVLT 30-second recall also correlated with functional connectivity of the left 
postcentral gyrus. Statistical maps were corrected at cluster level (20 voxels) across the whole brain. The scatter plots show the peak voxel correlations. Voxel-wise multiple comparisons 
are thresholded with 5% FDR correction. Statistical maps are superimposed on a rendering of the Montreal Neurological Institute template brain. The color scheme of each image is 
normalized to the peak voxel correlation with the respective neuropsychological score, p values on scatter plots indicate corrected p value for the peak voxel correlation. The corresponding 
rvalues and corrected p values at the 5% FDR threshold (rjthresh)) include: lexical fluency, r= 0.4825, p = 0.0093; category fluency, r = 0.4785. p = 0.01; digit span bacl(ward. r = 0.4795. 
p = 0.0098; CVLT30 s recall, r = 0.499. p = 0.0095; CVLT total score, r = 0.497, p = 0.0098. MCI = mild cognitive impairment, AD = Alzheimer's disease, PCA = posterior cortical atrophy. 



Amn/Dys = amnestic/dysexecutive, AD-Language = logopenic variant primaiy progressive aphasia, CVLT = 
threshold, r(max) = correlation coefficient of the peak voxel. 



California verbal learning test, rjthre 



: correlation coefficient at the 5% FDR 



describes residual variability due to measurement error and other un- 
measured factors. We specified a hierarchical prior for (i„: 

fit ~ N{li;, u\) 



MCMC chains from over-dispersed initial values. After 5000 iterations 
of burn-in, we obtained 5000 additional samples from each chain. Con- 
vergence was achieved, with all chains having Gelman Rubin statistics 
of 1 (Gelman and Rubin, 1992). We combined the three chains and 
thinned to obtain 5000 iterations with which to do inference. Mixing 
sufficed to yield > 2500 effectively independent samples (Plummer 
et al., 2006) from each chain. 



The voxel-specific effects ((i„) are centered around their test-level 
counterparts (fit), which are in turn centered around the domain-level 
effects (fif), which are in turn all centered around a common global (i. 
Domain-level effects correspond to the tests grouped as executive, vi- 
suospatial, and episodic memory categories, each reflecting a specific 
cognitive domain. The error terms are modeled as normal with mean 
zero and variance cr^e- The global p and all four as are given non- 
informative flat priors. The b, are modeled as N(0, x^). 

We fit the model with the Markov chain Monte Cario (MCMC) algo- 
rithm using the statistical programming language R (http://www.R- 
project.org ) to estimate the posterior distribution of each parameter. 
The 95% credible intervals represent the 2.5-97.5 percentiles of the pos- 
terior distribution. Additional multiple comparisons corrections are not 
required in the Bayesian hierarchical model since the shrinkage prior, 
which shifts estimates towards each other, results in 'partial pooling', 
a compromise between the overly simplistic global-level estimate and 
the overly noisy voxel-specific estimate. This method is well described 
in the statistical literature (Gelman et al., 2012). We started three 



3. Results 

3.3. AD-spectmm patients have reduced resting-state alpha activity and reduced 
pnctional connectivity ttiat correlates with higher disease severity 

Consistent with previous reports (de Haan et al., 2008; Jeong, 
2004), in a quantitative power analysis of MEG data, AD-spectrum 
patients showed significant slowing of alpha-band (8-12 Hz) ac- 
tivity compared to an age-matched control group (Fig. 1). To mea- 
sure functional connectivity, we computed imaginary coherence 
within the alpha band at each voxel (see Methods for details) 
depicting the connectivity of a given voxel to the rest of the brain. 
We found that functional connectivity is inversely related to CDR- 
SOB ( degree of disease severity) in AD spectrum. Higher CDR-SOB cor- 
related with reduced levels of imaginary coherence in clusters localized in 
the right posterior perisylvian region and left middle frontal cortex. The 
strongest cortelation was seen over the right posterior perisylvian region 
(Fig. 2 A and C) and the second strongest was over the left frontal cortex 
(Fig. 2B and D). These results indicate that neuronal synchronization of 
spatially distributed regions of the brain is significantly correlated with 
the global cognitive disability in AD spectrum. 
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3.2. Subgroups within AD spectrum have distinct patterns of cognitive 
dysfunction 

Neuropsychological measures of specific cognitive domains liave re- 
liably distinguished between distinct clinical variants of AD spectrum 
(Rabinovici et al., 2010). Consistent with these earlier reports, patients 
in the IvPPA, PCA and amnestic/dysexecutive subgroups in our study 
showed unique profiles of cognitive performance (Table 2). For exam- 
ple, PCA patients had better lexical fluency scores (e.g., D word genera- 
tion) than IvPPA patients (p < 0.05, one-way-ANOVA, Tukey post hoc), 
yet performed poorer in visual construction (e.g., Benson copy) than 
all other subgroups (p < 0.0001, one-way-ANOVA Tukey post hoc). 
Other measures that differed significantly among the clinical subgroups 
included digit span backward (worse in PCA and IvPPA than MCI), VOSP 
number location (worse in PCA than amnestic/dysexecutive variant and 
MCI), and Stroop color naming (worse in PCA than MCI). We predicted 
that resting-state functional connectivity deficits could account for im- 
pairments of specific cognitive functions in AD spectrum. Specifically 
we tested the hypothesis that poor cognitive performance would be 
correlated with low imaginary coherence measures of brain areas that 
are important for each task. 

3.3. Regional patterns of resting-state functional connectivity are correlated 
with focal cognitive impairments in AD spectrum 

With regards to executive function, we found that lexical fluency, 
categoiy fluency, digit span baclcward, short-delay verbal memory, 



and verbal learning were each significantly correlated with the func- 
tional connectivity of the left dorsolateral prefrontal cortex (Fig. 3A-E, 
Supplementary Table 2). Lower resting-state functional connectivity of 
this region predicted worse executive performance. IvPPA patients gen- 
erally scored low in lexical fluency and had low functional connectivity 
in the left dorsolateral prefrontal cortex, whereas PCA patients generally 
scored high on lexical fluency and had high functional connectivity in 
this region (Fig. 3A and, 3F triangles vs. diamonds). Interestingly, one 
subject met criteria for PCA yet had lower scores in lexical fluency 
than is typical for PCA; this individual showed relatively poor connectiv- 
ity of the left dorsolateral prefrontal cortex compared to other PCA sub- 
jects who performed better in lexical fluency (3F). This case illustrates 
associations between phenotypes and network connectivity that 
would not be predicted from clinical sub-classification alone. 

We next searched for brain regions where visuospatial function 
would correlate with MEGl functional connectivity in AD spectrum. Vi- 
sual construction ability (Benson copy) and location discrimination abil- 
ity (VOSP number location) correlated with the degree of resting-state 
connectivity of the left inferior parietal cortex (Fig. 4, Supplementary 
Table 2). Weak resting-state functional connectivity of the left inferior 
parietal cortex predicted poor visuospatial ability. PCA patients had 
low scores in Benson copy compared to IvPPA patients and likewise 
showed low functional connectivity over the left inferior parietal cortex 
in relation to IvPPA patients (Fig. 4C). 

Finally, we searched for brain regions where episodic memory func- 
tion would correlate with MEGI functional connectivity. Delayed verbal 
memory (CVLT 10-minute recall) defects correlated with connectivity 
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Fig. 4. Resting-state functional connectivity deficits in the left inferior parietal cortex correlate with impairments in visuospatial ability in AD spectrum. Resting-state functional connec- 
tivity of the left inferior parietal cortex correlated with performance on the spatial tasks of (A. C) visual construction (Benson copy), and (B, D) location discrimination (visual object and 
space perception (VOSP) number location). The scatter plots show the peak voxel correlations. Voxel-wise multiple comparisons are thresholded with 10% FDR correction. Statistical maps 
are superimposed on a rendering of the Montreal Neurological Institute template brain. The color scheme of each image is normalized to the peak voxel correlation with the respective 
neuropsychological score, p values on scatter plots indicate corrected p value for the peak voxel correlation. The corresponding r values and corrected p values at the 10% FDR threshold 
(r(thresh)): visual construction, r = 0.48, p = 0.0097; location discrimination, r = 0.51 55, p = 0.0099. MCI = mild cognitive impairment, AD = Alzheimer's disease, PO\ = posterior cortical 
atrophy, Amn/Dys = amnestic/dysexecutive, AD-Language = logopenic variant primary progressive aphasia, VOSP number location = number location task of Visual Object and Space 
Perception batteiy, r(thresh) = correlation coefficient at the 10% FDR threshold, r^max) = correlation coefficient of the peak voxel. 
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deficits in tlie left medial frontal association regions, whereas Benson 
recall defects correlated with connectivity deficits in the left lateral fron- 
tal association regions (Fig. 5A and C, Supplementary Table 2). The 
amnestic/dysexecutive subgroup showed the lowest scores in episodic 
memory and also the weakest functional connectivity of the left frontal 
association regions (Fig. 5B and D). 

Collectively, these results support our hypothesis that region- 
specific functional connectivity deficits are associated with specific cog- 
nitive defects in AD spectrum. Whole-brain analysis revealed no regions 
in which higher functional connectivity associated with worse perfor- 
mance on any of the neuropsychological tests assessed above. Addition- 
ally, we found no regions where functional connectivity correlated with 
language performance on the BNT. Given that AD-spectrum patients 
generally performed well on this task (Table 2), we were likely under- 
powered to detect the degree of effect size that was observed in the 
other cognitive domains. 

We performed a Bayesian validation analysis to 1 ) account for multi- 
ple testing and 2) ensure that the cluster of voxels around the peak 
voxel — and not just the voxel with the strongest correlation — 
contained a significant association with the neuropsychological score 
(see methods for details). The validation analysis quantified the effect of 
connectivity (imaginary coherence) of the voxels in each cluster on 
each of the cognitive test scores or CDR-SOB, and estimated the difference 
in each measure per unit change in imaginary coherence (Fig. 6). For ex- 
ample, the standard deviation of scores on the lexical fluency was 3.41. 



Therefore, the estimated effect of 0.25 indicates that, across voxels 
in that region, a 0.01 unit increase in imaginary coherence is asso- 
ciated with a 3.41 x 0.25 = 0.85 unit increase in the lexical fluency 
score. None of the confidence intervals (indicated by thin bars for 
each cognitive task) crossed zero, indicating that all estimated effects 
are statistically significant. This analysis revealed that after accounting 
for multiple comparisons and for the cluster of voxels in each region 
(as opposed to the peak voxel), there was still a significant relationship 
between functional connectivity and cognitive measures. The Bayesian 
analysis, depicted in Fig. 6, represents a unifying "meta-analysis" of 
our data across both spatial and cognitive domains. 

3.4. Regional atrophy patterns did not account for correlations between 
functional connectivity deficits and specific cognitive impairments in AD 
spectrum 

Neurodegenerative diseases cause loss of neurons and neuropil, 
which is accompanied by gray matter atrophy. The extent to which 
early cognitive and functional connectivity deficits directly relate to 
atrophy is unclear. Using a voxel based morphometric assessment we 
examined the gray matter atrophy patterns of AD. Consistent with pre- 
vious reports the PCA patients showed atrophy predominantly in 
temporo-parietal-occipital regions (Fig. 7A), while the IvPPA pa- 
tients showed atrophy predominantly over the left temporo-parietal 
cortex (Fig. 7C). The amnestic/dysexecutive subgroup showed mild- 
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Fig. 5. Resting-state functional connectivity deficits in frontal association regions correlate with impairments in episodic memory performance in AD spectrum. Resting-state functional 
connectivity of two distinct frontal association regions correlated with performance on the episodic memoiy tasks (A, B) CVLT delayed recall (verbal memoiy) and (CD) Benson recall 
(visual memory). CVLT delayed recall was related to a relatively medial region of the left frontal cortex (A), whereas Benson recall was related to a relatively more posterior and lateral 
region in the left frontal cortex (C). The scatter plots show the peak voxel correlations. Voxel-wise multiple comparisons are thresholded with 5% FDR correction. Statistical maps are 
superimposed on a rendering of the Montreal Neurological Institute template brain. The color scheme of each image is normalized to the peak voxel correlation with the respective 
neuropsychological score, p values on scatter plots indicate corrected p value for the peak voxel correlation. The corresponding r values and corrected p values at the 5% FDR threshold 
(r(thresh)): CVLT delayed recall, r = 0.498, p = 0.0096; Benson recall, r = 0.4795. p = 0.0098. MCI = mild cognitive impairment, AD = Alzheimer's disease, PCA = posterior cortical at- 



rophy, Amn/Dys = amnestic/dysexecutive, AD-Language = logopenic variant primary progressive aphasia, CVLT = 
5% FDR threshold, r(max) = correlation coefficient of the peak voxel. 
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to-moderate atrophy in the left hippocainpus and bilaterally in the 
teiTiporo-occipital and frontal cortices (Fig. 7B). The MCI subgroup 
(n = 4) was not significantly different froin controls. 

We searched for brain regions where gray matter atrophy would 
correlate with deficits in CDR-SOB or cognitive performance using an 
approach that was comparable to our functional connectivity analysis. 
We used the same tasks that were found to have significant correlations 
between cognitive impairments and functional connectivity defects. 
These included lexical fluency, category fluency, digit span backward, 
Benson copy, and VOSP nuinber location. We failed to identify any clus- 
ters that showed significant correlations between CDR-SOB or specific 
cognitive deficits and focal atrophy patterns (Supplementary Table 4). 
To more carefully scrutinize regions that were identified in our func- 
tional connectivity analysis, we pre-selected cortical regions that were 
localized in the MEGI analysis and searched for correlations between 
the gray matter volume estimates of these regions and CDR-SOB or re- 
lated cognitive function (Suppleinentary Fig. 1). This analysis also re- 
vealed no significant correlations between the volume estimates and 
the neuropsychological measures. Thus, in our cohort, the volume esti- 
mates obtained by MRl appeared to be less sensitive predictors of cogni- 
tive dysfunction than abnorinal neural synchronization patterns 
detected by MEGI. 

4. Discussion 

In this study, we showed direct correlations between specific cogni- 
tive iiTipairments and deficits in regional functional network connectiv- 
ity in AD spectrum. Several unique features of our investigation enabled 
these novel observations: 1 ) our methodology was unbiased and did not 
preselect neural networks to interrogate, 2) we performed the MEG 
analysis in source space (brain space) rather than sensor space enabling 
us to examine brain-region specific connectivity patterns, 3) we directiy 
measured electrodynamic changes in alpha activity rather than low- 
frequency hemodynamic changes that are measured with fMRI, and 4) 




Fig. 7. Voxel based morphometry (VBM)-derived atrophy patterns for different clinical 
variants of AD. VBM atrophy maps based on comparisons with age-matched healthy con- 
trols are shown for the three clinical variants of Alzheimer's disease: (A) posterior cortical 
atrophy (PCA), (B) amnestic/dysexecutive subgroup (Amn/Dys). and (C) logopenic vari- 
ant primary progressive aphasia (AD-Language). Regions of gray matter atrophy are 
shown on the 3-dimensional rendering of the Montreal Neurological Institute (MNl) stan- 
dard template brain. Results for PCA and AD-Language were corrected for family-wise 
error (p < 0.05) and the results for the amnestic/dysexecutive subgroup were thesholded 
for uncorrected p < 0.001 . MNl coordinates and corresponding t values are provided in 
Supplementary Table 3. n = 7 PCA, n = 7 amnestic/dysexecutive, and n = 4 AD-Language. 



our study population comprised three clinical subtypes of AD as well 
as MCI, which enhanced the sample heterogeneity and increased the 
likelihood of identifying distinct patterns of network dysfunction. The 
degree of gray matter atrophy in our cohort did not correlate with cog- 
nitive dysfunction. Although a direct comparison of the MEGI-derived 
functional connectivity index and the MRI-derived gray inatter volume 
index is beyond the scope of this paper, stronger correlations between 
MEGI and cognitive perforinance suggest that abnormalities in 
resting-state functional connectivity are more sensitive indicators of 
early cognitive dysfunction than atrophy patterns. Our results deinon- 
strate that region-specific neural dysconnectivity patterns are a poten- 
tially useful and precise tool to map and follow defective networks in 
the early stages of AD. 

These findings complement and extend upon previous electrophys- 
iological studies of AD. Tools that allow mapping of distributed function- 
al networks such as EEC and MEG have demonstrated distinct patterns 
of abnormal connectivity in AD spectmm (Brenner et al., 1986; de Haan 
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et al., 2008; de Haan et al., 2012b; Duffy et al., 1984; Martin-Loeches 
et al., 1991 ). Previous reports revealed predominantiy left hemispheric 
decreased global alpha and beta-band connectivity and increased parie- 
tal theta connectivity in AD (de Haan et al., 2012a; de Haan et al., 2008; 
de Haan et al., 201 2b; Osipova et al., 2005; Stam et al., 2009; Stam et al., 
2006). The current results demonstrate that focal patterns of alpha 
dysconnectivity directly relate to domain-specific cognitive impairments 
in AD spectmm. 

Our unbiased MEGl approach identified brain regions that have been 
implicated to relate to behavior in previous functional and anatomical 
studies. The dorsolateral prefrontal cortex forms part of a network 
that is involved in executive performance, attention, working memory, 
and learning (Dosenbach et al., 2007; Seeley et al., 2007; Vincent et al., 
2008; Yeo et al., 2011). Likewise, parietal cortical regions have been 
linked to visuospatial function (Cabeza and Nyberg, 2000). Consistent 
with these findings, our results linked connectivity deficits of the left 
dorsolateral prefrontal cortex to executive dysfunction and connectivity 
deficits of the left inferior parietal cortex to visuospatial dysfunction. We 
also found that reduced functional connectivity of left frontal associa- 
tion regions are correlated with episodic memoiy deficits. This finding 
is consistent with previous experiments showing that frontal lobe dam- 
age impairs both encoding and retrieval of memory (Cabeza and 
Nyberg, 2000; Goldman-Rakic, 1987, 1988; Creshberg and Shimamura, 
1995; Janowsl<y et al., 1989). It is noteworthy, though, that the medial 
temporal lobes including the hippocampi, which play an important role 
in memory, were poorly measured in the present analysis because the 
current MEGl technology does not allow reliable interrogation of deeper 
brain structures. Algorithms designed to explore signals from deeper 
brain structures would resolve this constraint in the future. 

We found overlapping brain regions of functional connectivity that 
correlated with both category and lexical fluency. This observation is 
contrary to the widely accepted notion that these two tasks are related 
to different anatomical regions. Previous studies have demonstrated 
that lexical fluency primarily involves the frontal cortex, whereas cate- 
gory fluency primarily involves the temporal cortex. Some of the most 
compelling evidence for this distinction comes from investigations of 
normal aging populations, patients with semantic dementia caused by 
frontotemporal lobar degeneration, and stroke patients (Baldo et al., 
2006; Grogan et al., 2009; Libon et al., 2009). However, other studies 
have reported that both lexical and categoiy fluency are similarly im- 
paired in patients with frontal lesions (Baldo and Shimamura, 1998; 
Schwartz and Baldo, 2001). Our results support the hypothesis that 
patients with frontal lobe pathology have a reduced ability to make 
strategic and effective searches through memory, independent of 
whether the search is phonemically or semantically driven (Hirshorn 
and Thompson-Schill, 2006; Troyer et al., 1 998). The mechanism under- 
lying poor categoiy fluency in AD could result from the defective net- 
work connectivity of frontal cortices, rather than from an inherent 
defect of the semantic association regions of the temporal cortex as oc- 
curs in semantic dementia. 

Future studies are needed to further characterize the disparate net- 
works involved in AD and to address some of the limitations of this 
study as follows. First, the sample numbers of specific AD-spectrum var- 
iants were relatively small. Hence, group contrasts between clinical var- 
iants could not be performed. Second, we focused only on the alpha 
frequency range because there are currently no reliable or reproducible 
methods to reconstruct source-space functional connectivity in other 
frequency bands, such as beta or gamma. Third, the patients seen at 
our specialized center included mostly early-onset disease. Although 
this approach enabled us to identify specific AD variants, it also limited 
generalization to the more commonly encountered late-onset AD popu- 
lation. Fourth, while the biomarkers strongly support that our cohort 
consists of AD pathology, we do not have autopsy confirmation. Future 
investigations recruiting more numbers, extending to the whole 
spectrum of AD, with follow-up autopsy confirmations, and with addi- 
tional analyses of other frequency bands, will further advance our 



understanding of network dysfunctions in AD spectrum. Additionally, 
an a priori region-of-interest analysis in a new cohort using the regions 
and correlated tasks that we identified, would be an ideal model for a 
follow-up study. 

In conclusion, the current study details connectivity deficits in AD 
spectrum by measuring the direct electrical activity of the brain with 
MEGl. The strength of regional network connectivity is significantly cor- 
related with the performance in specific cognitive domains. MEGl func- 
tional connectivity provides a new approach to study dysfunctional 
networks in AD and other neurodegenerative dementias with high spa- 
tial and temporal resolution. 
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